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Abstract. We extend the results obtained in [T2] about a class of Lagrangian sys- 
tems which admit alternative kinetic energy metrics to second-order mechanical sys- 
tems with explicit time-dependence. The main results are that a time-dependent al- 
ternative metric will have constant eigenvalues, and will give rise to a time-dependent 
coordinate transformation which partially decouples the system. 



1 Introduction 

In a recent contribution [12] two of us examined the following particular case of the 
inverse problem of Lagrangian mechanics: suppose that an autonomous second-order 
system admits a Lagrangian description such that the Hessian of the Lagrangian L is the 
unit matrix; what are the conditions for the existence of an alternative kinetic energy 
metric for the same system, or expressed differently, an alternative multiplier for the 
inverse problem, which depends on the position variables only. Two observations we 
made are worth mentioning right away: (i) the starting assumption implies that the 
forces under consideration necessarily are of (generalized) electromagnetic type; (ii) if 
a second multiplier exists, it necessarily will have to be constant. We obtained several 
interesting results for such systems, one of them being that they necessarily decouple in 
coordinates which diagonalize the alternative metric. 

It turns out to be quite interesting to enlarge the scope of applicability of the preceding 
study by allowing explicit time-dependence into the picture. This is the subject of the 
present paper. We will see indeed that such an extension leads to non-trivial complica- 
tions, such as the fact that the determination of admissible scalar and vector potentials 
is no longer a separate issue. Nevertheless, the conclusion about partial decoupling of 
the system in appropriate coordinates will be shown to remain intact. 
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The literature about the inverse problem of Lagrangian dynamics is vast. For a small 
sample of different techniques and mathematical tools which have been used, see [10\ [21 
El El EE] and the list of references in the review paper [SJ. A number of these techniques 
require quite advanced methods of differential geometry. Note that allowing explicit 
time-dependence in the framework is not always a straightforward matter; it involves, 
for example, passing from the vector bundle structure of a tangent bundle to the affine 
context of a first jet bundle. For specific aspects of time-dependence in this context, see 
for example [U [13] . The very concrete problem we are addressing here is much more 
computational in nature, so that it will be sufficient to rely on analytical expressions of 
the conditions for the existence of a Lagrangian, such as the ones developed in |10| . 



The second-order systems (sodes for short) we start from are general non-autonomous 
differential equations in normal form, say 

x a = F a (t,x,x), a = l,...,n. (1) 

Geometrically they are the analytical expression of a second-order equation field T living 
on the first jet bundle J l ir of a bundle 7r : E — > M, so 

r -£ + **5? + **<«■ **>Si- 

The local formulation of the general inverse problem is the question for existence of a 
(non-singular) multiplier matrix g ab (t, x, x), such that 

■ b ^ d ( dL\ dL 



g ab (x b — F ) = — 



dt \ dx a J dx a ' 

for some Lagrangian function L(t,x,x). The necessary and sufficient conditions for the 
existence of L are generally referred to as the Helmholtz conditions, and when expressed 
in terms of conditions on the multiplier matrix, they read (see 

9ab = 9ba, F(g ab ) = g ac T c b + g bc T c a , g ac <& c b = g bc $ c a , 

Here 

ldF a 



dg a b dg a 



dx c dx b 



-pa 

b '~ 2 8x b 

geometrically represent the components of the canonical non-linear connection associated 
to r, so they are non-tensorial objects. On the other hand, 

dF a 

n :=-^r-rgr«-r(r£) 

are the components of a type (1, 1) tensor field along the projection ttq : J 1 ^ — >■ E called 
the Jacobi endomorphism. For a systematic development of the geometrical calculus 
along 7Tq, the interested reader can consult [13]. Further worth mentioning for later 
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reference is that an integrability study of the conditions for the g ab leads to more algebraic 
restrictions, the first one of which, called the curvature condition, reads 

QarKc + QbrKa + 9crKb = °- ( 3 ) 

Here, 

R a bc :=H c (T a b )-H b (T a c ), where H c = JL - r^A . 

are components of the curvature tensor of the non-linear connection referred to before. 

The concrete question we shall address in this paper is the characterization of (generally 
time-dependent) second-order systems which have a Lagrangian representation with a 
Euclidean kinetic energy term, and which admit at least one alternative Lagrangian 
representation for which the multiplier matrix does not depend on the velocities but 
may explicitly depend on time. The main conditions to be satisfied by such systems 
are developed in section 2; one of the conclusions is that the alternative multiplier g 
can in fact be a function of time only. In section 3, it is shown that this g will have 
constant eigenvalues and that a linear coordinate transformation which diagonalizes g 
will decouple the system into lower dimensional Lagrangian systems. Explicit examples 
which illustrate the various aspects of the theory are presented in section 4. In section 5, 
we are led to an additional characterization of the systems under consideration, in the 
particular case that the connection coefficients are functions of time only. 



2 Lagrangians of electromagnetic type and alternative mul- 
tipliers 

Suppose r is a SODE field on J l ir with the property that, in some coordinate system, 
the unit matrix: g ab = 5 ab is a multiplier matrix for the inverse problem. In fact, for 
simplicity, we shall take tt to be the trivial bundle ir : R x IR n — > R. Since the unit 
matrix is the Hessian of the function ^ ^(i a ) 2 , our assumption means that there exists 
a Lagrangian L for the given SODE, which differs from that quadratic function at most by 
a linear function in the velocities, say L = ^ ^2(x a ) 2 + A c (t, x)x c — V(x, t). This implies 
that the 'forces' will be of the form (A a = 5 ab A b ) 

pa= (dA^_dA«\ ±b _dV_dA^ 

\ dx a dx b dx a dt ' K ' 



va _ _ l dF a _ 1 (dA a dA b 



It further follows that the connection coefficients are skew-symmetric: 

(5) 

while the components of the Jacobi endomorphism become 

^a_\( d 2 A a 8 2 A C d 2 A b \ . b ^ d 2 V 1 ( d 2 A a 8 2 A C \ 
c ~ 2 \dx c dx b + dx a dx b ~ dx a dx c ) X ~ b c+ dx a dx c + 2 \dtdx c + dtdx a J ^' 
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and are seen to be symmetric. 



We want to investigate under what conditions there exists a second multiplier for such 
systems, which does not depend on the velocities. If g does not depend on the velocities 
it follows from the requirement 



F{g a b) 



7« b 



+ 



Jab ■ c 
-X 



9ac^b + 9bcFa-> 



dt dx c 

that g must actually be a function of t only for which it must further hold that 



9ab(t)=g ac (t)T$ + g bc (t)T c a . 



(7) 



Notice that this implies 



o. 



(8) 



which is a somewhat weaker skew-symmetry property than in the autonomous case. The 
remaining Helmholtz condition g a c^1 = 9bc& a 1S a linear expression in the x a as well and 
thus also splits into two conditions. From the terms which are linear in the velocities 
it can be shown that the weaker condition (jHJ) is sufficient to imply that (as in the 
autonomous case) we still must have 



djar dj rb dj ba 

dx b dx a dx r 



0: 



(9) 



with := gac^t- This may surprise at first sight, but it is fully consistent with the 
curvature condition (|3|) for this case. The remaining terms in the ^-condition, which in 
the autonomous case produced conditions on admissible scalar potentials V only, now 
give rise to the following more complicated requirements, involving both the scalar and 
vector potential: 



9a 



d 2 V t 1 f d 2 A c d 2 A b 
dx c dx b + ^ \dtdx h + dtdx' 



■pc-pr 



9bc 



d 2 v , 

h - 

dx c dx a 2 



Q2 A c Q2 A a 
+ 



dtdx a dtdx 1 



(10) 



Summarizing, we come to the following conclusion. 

Proposition 1. Assume that the SODE T on R 2n+1 admits a coordinate representation 
for which 5 a b is a multiplier matrix for the inverse problem of Lagrangian mechanics. 
Then the right-hand sides of the equations are of the form for some functions A a and 
V . Any admitted alternative multiplier g a b, not a multiple of the identity and independent 
of the velocities, must be a function of t only, must satisfy the differential conditions $^), 
the curvature type conditions ^) with respect to the functions A a , and the conditions \10\) 
with respect to the functions V and A a . □ 
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In the next section, we shall see that these conditions imply existence of a linear change 
of coordinates which partially decouples the given equations. 

An interesting special case which complements the results of [12] is when the system 
with given forces (J4j) does depend explicitly on time, but we restrict the search for an 
alternative multiplier to a constant g. It then follows from ([7|) that the matrix gT must 
be skew-symmetric. As a result the terms with two factors T in the condition (jlOp cancel 
out, but compared to the autonomous case in [12] there remains a coupling between the 
scalar and vector potential. From the fact that also gdT/dt will be skew-symmetric then, 
the remaining condition (|1(J|) can be further manipulated to reduce to 

d fdV 8A C \ _ d fdV dA a \ 
^ ac dx b \dx c dt J ^ bc dx c \dx a dt ) 



3 Diagonalization of g and decoupling 

We observe that the differential conditions (JT]) for g have the form of a Lax equation 
because the matrix is skew-symmetric. It follows that the eigenvalues of g are constant. 
The point now is that the diagonalization of g can be achieved by a linear coordinate 
change and that the connection coefficients, although not behaving tensorially under such 
operation, transform exactly in the way which is needed to contribute to decoupling of 
the system. 

Let -P(t) be an orthogonal matrix which diagonalizes the symmetric matrix g(t) and put 
g = P T gP, so that g is diagonal. Using an obvious matrix notation, the equation ([7]) 
can be written as 

g = gT - r#. 

[It will always be clear from the context whether the notation T refers to the SODE vector 
field or to the matrix of connection coefficients.] Multiplying on the left by P T and on 
the right by P, the equation transforms to 

g = g(T + P T P) - (f + P T P)g, (11) 

where f = P T TP. Notice that T is still skew-symmetric and the same is true for 
r := r + P T P. It trivially follows from the fact that g is constant that the commutator 
in the right-hand side is zero. 

Now consider the linear, time-dependent coordinate transformation x = P(t)y with 
inverse y = P(t) T x. Both g a ^ and are actually components of tensor fields along 
7Tq : J 1 it — > E of the form (see |13j ) 

g = g ab e a ®6 b , $ = $ a b 9 b with 6 a = dx a - x a dt. 

ox a 

Hence they behave well under time-dependent coordinate transformations, meaning that 
their components in the new variables under consideration will be g and <l = P T QP. 
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Less obvious is what happens to the non-tensorial matrix of connection coefficients. It 
is a straightforward computation, however, to verify that the transformed second-order 
differential equations for the y- variables contain the following linear part in the velocities: 

y = -2{P T TP - P T P)y + .... 

Thus, the connection coefficients T% of the transformed equations are precisely the com- 
ponents of the matrix V introduced above. Since T commutes with the diagonal matrix 
g, it has a block diagonal structure, the dimension of each block being determined by the 
multiplicity of an eigenvalue of g. It follows that, if A a and Bp refer to coordinates of 
different blocks, the transformed forces F will have the property dF Aa /dy B ^ = 0. But 
we also have g& = <i>g, meaning that <I> has the same block diagonal structure, and from 
the coordinate expression of (see ([6]), but to be read in the new variables), it then 
further follows that also dF Aa /dy B i 3 = 0. Hence, the transformed equations decouple 
into independent subsystems, one for each eigenspace of g. 

Proposition 2. Assume that a given SODE T satisfies the conditions of proposition 1. 
Then the linear time- dependent coordinate transformation which diagonalizes the second 
multiplier g a b(t) will decouple the system into a number of separate SODEs, one for each 
distinct eigenvalue of g, and these eigenvalues are constants. □ 

As we explained in [12] such a decoupling property was already established, for au- 
tonomous systems, in work of Ferrario et al [5j, but to the best of our knowledge, it 
was not observed before that allowing time-dependence in the picture does not destroy 
it. Observe also that the transformed system (and hence every independent subsystem) 
still has the property of admitting the unit matrix as multiplier, because P T I n P = I n . 
As a result, a reverse engineering technique to construct all systems under consideration 
goes as follows: start from a number of say k independent SODEs with forces of type (HJ); 
multiply each of them with a different constant factor thus creating a diagonal multi- 
plier g for the union of the k systems; then apply an arbitrary (possibly time-dependent) 
linear orthogonal transformation to the overall system; the resulting coupled system will 
be of the type characterized in Proposition 1. Note that each choice of a single equation 
at the start of this process, corresponding to an eigenvalue with multiplicity one of the 
g under construction, will only contain a scalar potential. Since the group of orthogonal 
transformations is finitely generated, one can imagine that such a reverse engineering 
technique may be helpful in an attempt to classify all systems in the category under con- 
sideration. But there would of course, in such a classification, be room for an unlimited 
number of ways of putting together decoupled systems in the first place. 

The more interesting approach for identifying suitable examples, therefore, is to start 
from an arbitrary SODE with forces (JH) and to try to solve the conditions of Proposition 1 
for the identification of admissible scalar and vector potentials and the construction of 
an associated alternative (possibly time-dependent) kinetic energy function. This is the 
line of approach we will adopt as much as possible in the next section. 



6 



4 Examples 



Consider a system x a = F a {x, x,t), where the F a are of the form ([4]). In the autonomous 
case we were able to obtain relevant information about the form of admissible vector 
potentials and the alternative multiplier g by imposing first that the matrix 7 := gT 
had to be skew and satisfy the curvature type condition Q. The remaining Helmholtz 
condition g& = (g&) T then resulted essentially in equations for the admissible scalar 
potential V. The situation is more complicated in the time-dependent case now, for 
the following three reasons: firstly, 7 need not be skew, instead, in view of dH), we 
have the weaker condition that 7 — j T must be a function of time only; secondly, the 
information about g which can be gathered in the process of analysing 7 will be far less 
conclusive, because g still has to satisfy the differential equation ([7]) in the end; thirdly, 
the remaining algebraic condition (jlOp will in general not reduce to a separate issue 
concerning admissible functions V. As we shall show below, however, we can still get a 
long way in constructing examples for n = 2 and n = 3 by staying close to the procedure 
we followed in the autonomous case. 

Taking n = 2, and denoting the elements of the multiplier g for shorthand by at = 
9n, 02 = 522, b = gu = #21 1 the fact that 7 — 7 T can only depend on t implies that 
bY\ and (ai — 02)1^2 must be functions of time only. It follows that Y\ cannot depend 
on the ^-coordinates either (which we write with lower indices here); putting T\ = <r(t), 
the vector potential must satisfy 

8x2 dxi ' 
and is determined, up to a gradient, by 

Ax{t, x) = a(t) x 2 , A 2 (t, x) = -a{t) x x . (12) 

The freedom of adding a gradient is irrelevant here since this adds to the Lagrangian L 
we started from merely a total-time derivative term, plus a term which can be absorbed 
into the as yet undetermined scalar potential V. Expressed differently, we can assume 
without loss of generality that in dimension n = 2 the system we start from can be 
written in the form 

dV 

x t = -2a(t)x 2 - ox 2 - 3 — , 

ox 1 

dV 

'±2 = 2a(t)xi + vxi - t; — • 

0x2 

Quite remarkably, if we now move to the remaining algebraic condition (jlOp . all terms 
not involving V cancel out and we end up with the following condition for admissible 
functions V: 

'd 2 V d 2 V\ , d 2 V 

(02 ~ 01) „ , . (13) 



dx\ dx\ ) dx\dx2 
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By way of example we proceed to determine admissible (time-dependent) quadratic po- 
tentials. Putting 

V = \k(t)x\ + l(t)x 1 x 2 + \m{t)x\, 

the condition becomes 

b(m — k) = (a 2 — a\) I 
and we will only discuss the generic case here that I 7^ and m 7^ k. 



Putting l(t) = (m(t) — k(t))p(t) the multiplier must have the form 

9 = 



at (a 2 -ai)p 
(a 2 -a{)p a 2 

and the corresponding admissible potentials read 

V = \kx\ + (m — k)px\x 2 + \mx\. (14) 

But now g still must satisfy the differential equations ([7]). Putting a 2 = C — a\, where 
C represents the constant trace of g, these equations reduce to: 

di = 2(2ai — C)pa 
p = _(i + 4/0 2 )a. 

The case 2a\ — C = has to be excluded since g must be kept different from a multiple of 
the identity. It is convenient to introduce a new time-scale r := J a(t)dt; the equations 
then are easy to integrate, giving 

p(r) = |tan(-2r + ^), 2a 1 {r) - C = B cos{-2t + A), 

where A and B are constants. 

In conclusion, with k(t), l(t), a{t) arbitrary functions of time, a(t) = —2f t a(s)ds + A 
and p(t) = 2 tana, the admissible scalar and vector potentials are given by (TT4"|) and (fT2l) 
respectively, and the alternative multiplier is 

\{C + B cos a) —TjBsma 

I 

2 



1 Bsina \{C — B cos a) 



The eigenvalues of g are Ai = ^(C — B), X 2 = |(C + B) and corresponding eigenvectors 
read (sin a, cos a + 1), (sin a, cos a — 1). They have to be scaled to norm 1 to produce 
the columns of an orthogonal matrix which diagonalizes g. For testing the statement 
about decoupling in Proposition 2, however, one can relax the requirement that P be 
orthogonal. Indeed, if we look at the related coordinate transformation in the form 
y = P T x whereby the eigenvectors now appear as rows, it is clear that an overall factor 
in the defining relation of the y-variables cannot be an obstruction to decoupling. One 
can verify that the transformation 

yi = sin a x\ + (cos a + \}x 2 
y 2 = sinaxi + (cos a — l)x 2 , 
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effectively takes care of that. 

Let's move now to systems with three degrees of freedom, for which the problem is more 
intricate because the curvature condition ([9]) comes into play. To better illustrate the 
effect of allowing time-dependence we try to stay close here to the analysis and notations 
of the n = 3 example in [12]. As indicated before, the matrix 'y = gT need not be skew- 
symmetric now, but its symmetric part, according to ([8]) can depend on time only. In 
particular, from the diagonal elements of this symmetric part, the following expressions 
(which had to be zero in the autonomous case) now must be certain functions of time 
(which we will not specify at the moment): 

917^2 + 9X3^1, 512^-523^, 513^+523^, 

and we can view this as a linear system for the T^(t, x) (with time-dependent coefficients) 
which has zero determinant. Except for special cases where the rank of the coefficient 
matrix is lower than two (and which we will not discuss), solving this algebraically for 
the F% will give rise to relations of the type 

T\ = a{t)T\ + ~a{t), V 2 = b{t)T\ +b(t), (15) 

for certain functions a,b,a,b. The sort of privileged role which Y\(t,x) plays in these 
expressions is just a matter of renumbering the variables, if necessary. Substituting this 
structural information back into the linear equations from which it was obtained, we see 
that the product of T\(t,x) with any of the factors 

512 + ^13, 9i2-bg 2 3, ag 13 + bg 2 3, 

must be some function of time only. There are then two cases to be considered: either 
r3, also depends on time only, and then the same is true for Tg and r|; or we insist that 
1^ should depend on at least one of the x-coordinates, and then the above mentioned 
factors necessarily must be zero. We investigate the second option first. 

With Tg and T| as specified above, the off diagonal elements of the symmetric part of 7 
give us three more algebraic conditions, which are affine linear in T^. Since we assume 
x-dependence in T^, such relations can hold true only if the coefficients of vanish 
which yields 

9u - 922 = bgi 3 + ag 2 3, a(gn - g 33 ) = g 23 - bg 12 , b(g 22 - g 33 ) = -ag 12 - g 13 . 

We now have sufficient elements to get an idea of the possible structure of the multiplier 
g. We must have g\ 2 = —ag±3 = bg 23 and to make ag\ 3 + bg 23 zero as well, we can put 
513 = —pb,g 23 = pa, which implies that g\ 2 = pab. Putting further g\\ = pc (which 
defines the function c(i)), the last three conditions fix g 22 and 533. The result is that the 
multiplier g will have the structure: 

(c ab —b\ 

ab c + b 2 - a 2 a , (16) 

-b a c + b 2 -1 / 
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where a, b, c and p are all functions of time which will have to be determined from the 
differential equations which g must satisfy. It is instructive at this stage to compare with 
the situation in the autonomous case (see |12j). It is not true here that if g is a multiplier 
matrix, then p(t)g is one too, so that the overall factor cannot be disregarded right away. 

In dimension three the curvature type condition contains one requirement and with the 
information gathered so far, it can be seen to reduce to 



a 7-^ + 6tt^ = 0, 



dx 3 8x2 dx\ 

which is formally the same as in the autonomous case. The fact that the coefficients 
depend on time now is not an obstruction to take over suitably adapted results we 
know from the autonomous situation. The general solution of the above linear partial 
differential equation is an arbitrary function = f(t,u,v) where 

u(t,x) = x\ — b(t) x$, v(t, x) = X2 + a(t) 23. (17) 

This means that admissible vector potentials must satisfy 

: 2f(t,u,v), 

-■ 2a(t) f(t,u,v) + 2a(t), 

2b(t)f(t,u,v) + 2b(t). 

The right-hand sides are the components of a closed 2-form indeed (with t treated as 
parameter), and a particular solution for the Ai then is given by 

A\ = 2v sf(t, su, sv)ds + d(t)x 3 , 
Jo 

A 2 = -2u sf(t,su,sv)ds + b(t)x 3 , (18) 
Jo 

A3 = —2(au + bv) / sf(t,su,sv)ds — a(t)xi — b(t)x2- 



dA 1 


8A 2 


dx 2 


dxi 


dAi 


dA 3 


dx 3 


dx\ 


8A 2 


dA 3 


dx 3 


dx 2 



JO 

The general solution differs from this particular one by a gradient, but as argued in the 
preceding example, this gradient term can be incorporated into the as yet undetermined 
scalar potential, up to a total time derivative. We have obtained valuable information 
about the structure of g and the vector potential now, but it is clear that the freedom 
in selecting functions a, b, c and p will be severely restricted if we now go back to the 
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differential conditions (|7|). They read: 



-^(pc) = 2pab, 

±-(p(c + b 2 -a 2 )) = -2pba, 
at 

±(p(c + b 2 )) = 2p(ba-ab), 

— (pab) = p{bb — aa), 
-^(pa) = p(b(l - a 2 ) + aab), 

— (-nh\ = n(nn - h 2 ^ 

dt 



.(-pb) = p(a{l-b 2 ) + bab). 



It follows from the second and third equation that we must have 

j t {pa 2 ) = 2p(-ab + 2ba). 

Combining this with information coming from the fifth and sixth equation, one can derive 
the condition 

ab — ba = 0, 

and subsequently that 

a = b and b = —a, 

which in turn implies that a 2 + b 2 is constant, and from the third equation also that 
p(c + b 2 ) is constant. Finally, going back to the fifth and sixth equation, which reduce to 

|M=A |(P6) = -PS, 

and using the equations for a and b just obtained, it follows that actually also p must be 
constant and can therefore, from now on, without loss of generality be put equal to 1. 
With these restrictions, all six equations about g are now satisfied. We conclude that the 
freedom in constructing alternative multipliers consists of two arbitrary constants and 
one arbitrary function of time. Take a(t) to be arbitrary for example, and choose two 
constants c\ and C2, then define b by a 2 + b 2 = c\ and c by c + b 2 = 02- This fixes g, and 
the corresponding admissible vector potentials are now also determined by the relations 
(fT8|h where f(t,u,v) is an arbitrary function of the indicated variables and d(t) = —b, 
b(t) = a. 

It remains to find the form of the admissible scalar potentials for which purpose we 
will avoid the conditions (I10p which are quite complicated to handle. Observe first 
that g (as in (fT6|) with p = 1) has a double eigenvalue A = c + b 2 = C2 and a single 
eigenvalue p, = c — a 2 — 1 = C2 — c\ — 1. From a different perspective, the solution 
we obtained for admissible vector potentials more or less forces a kind of privileged 
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coordinate transformation upon us. Indeed, extending (|17p . consider the coordinate 
transformation (similar to the autonomous case in |12j ) 



u 



v 



z 



xi - b(t)x 3 , 
x 2 + a(t)x 3 , 
x 3 + b(t)xi - a(t)x2- 



(19) 



The inverse transformation reads 



xi = k^ 1 ((a 2 + 1) u + abv + bz) 
X2 = k~ l (abu + (b 2 + 1) v — a z) 
x 3 = k~ x (—bu + av + z), 



(20) 



where k = a? + b 2 + 1. If P(t) denotes the coefficient matrix in these inverse relations, 
it is easy to verify that its first two columns are eigenvectors of g corresponding to the 
eigenvalue A, and the third column is an eigenvector for the eigenvalue fx. Although this 
P is not an orthogonal matrix, it is still true that P _1 gP diagonalizes g. Therefore, 
Proposition 2 guarantees that the equations of motion will decouple into a 2-dimensional 
system for u and v and a separate second-order equation for z. Reversing the argument, 
this implies that admissible scalar potentials for our original system necessarily must be 
of the form 



U and Z being arbitrary functions of the indicated variables. This concludes the analysis 
for the case that in the relations (fT5j) was required to depend explicitly on at least 
one of the x-variables. 

What about the case that Tr, depends on time only, say = cr{t)l Using the represen- 
tation (|15p . the matrix of connection coefficients then has the form 



As argued before, we can fix the admissible vector potentials, without loss of generality, 
to be 



since the addition of gradient terms can be absorbed in the scalar potential up to a total 
time-derivative. The curvature condition Q is automatically satisfied here but, unlike 
the situation we encountered in our discussion of examples for n = 2, the remaining 
conditions (|10|) do not reduce to equations involving the function V only and in fact 
turn out to be quite complicated. A different line of approach can be adopted for this 
situation, which in fact applies to systems with arbitrary dimension and is explained in 
the next section. 






(22) 




(23) 
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5 The case of connection coefficients depending on time 
only 



The inspiration comes from the discussion on linear systems in [TT]: it is a well known 
practice for the integration of systems of linear second-order equations to eliminate the 
terms linear in the velocities by an appropriate coordinate transformation. This idea 
applies just as well for non-linear equations as long as they have linear dependence on 
the velocities. If the T% are functions of time only, the equations we start from, written 
with matrix notations, are 

dV 

x + 2T(t)x + —(t,x) + T(t)x = 0. (24) 
ox 

If U(t) is the solution of the matrix differential equation 

U + TU = 0, U(0) = I n , (25) 

then g(t) = (U~ l ) T SU^ 1 is the solution of the equation g = gT + {gT) T with initial value 
g(0) = S, and the coordinate transformation x = U(t)y will eliminate the (linear) velocity 
dependence in the forces. But remember that T is skew-symmetric in our case, so that 
(I25p implies that U T — U T T = 0, from which it subsequently follows that U T U +U T U = 0. 
Hence U T U is constant and since f7(0) = I n , U(t) will actually be an orthogonal matrix 
with the identity matrix as initial value. It is further straightforward to verify that the 
transformed differential equations for the y-variables will read: 

BV 

y - {U T T 2 U)y + — = 0, with V(t,y) = V{t, U(t)y), 

or 

y + ^ = 0, with W(t,y) = V -\y T {U T T 2 U)y. (26) 

Observe that the y-equations still have the unit matrix as multiplier since the forces are 
of potential type. More interestingly is the fact that the alternative multiplier g(t) for 
the original equations translates to the existence of a constant alternative multiplier S for 
the y-equations. What's more in making this transition we have eliminated the problem 
that the remaining conditions (|10p had the complicating coupling between scalar and 
vector potential, because they now simply read 

d 2 W d 2 W 
^ ac QycQyb ^ bc dy c dy a 



We summarize our results in the following statement. 

Proposition 3. Assume that the SODE T on R 2n+1 admits a coordinate representation 
for which 5 a b is a multiplier matrix for the inverse problem of Lagrangian mechanics 
and that the connection coefficients are functions of time only. Then, without loss of 
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generality, the vector potential can be taken to be of the form A a (t,x) = T%(t)x , and for 
there to exist an alternative multiplier, it suffices to construct a non-singular, constant 
symmetric matrix S (not a multiple of the identity) and a function W(t,y), satisfying 
the requirements (27\ ). IfU(t) is the solution of the matrix differential equation [25\), the 
admissible scalar potentials for the given system will be of the form 

V(t, x) = W(t, U T {t)x) + \x T T 2 (t)x, (28) 

and the corresponding alternative multiplier is given by g(t) = U(t)SU T (t). □ 

It should be emphasized that the orthogonal transformation with the matrix U in this 
picture is different from the one with P used before: P is dictated by the diagonalization 
of the multiplier g, whereas U is dictated by the matrix of connection coefficients T(t). 
In other words there is no reason why the constant matrix S in Proposition 3 (which is 
the initial value of g) would be diagonal and therefore, the transformed equations (I26D 
are different from the decoupled equations which Proposition 2 guarantees to exist. One 
can easily combine the results of both propositions: once an admissible system would 
have been reduced to the form (|26p and we accordingly know the initial value S of the 
alternative multiplier g, it suffices to diagonalize S to find a further transformation which 
will decouple the system. Conversely, starting from a number of decoupled systems, each 
of the form zu\ = —dW^/dz for some scalar potential W^, multiply each subsystem 
with an appropriate constant Aj and then apply an arbitrary constant, orthogonal trans- 
formation to the direct sum total system; the result will be a potential system which has 
the property (|27|) . with W = expressed in the new y-variables. Next, select an 

arbitrary orthogonal matrix U(t) which has the unit matrix as initial value and define 
T(t) by T := —UU T . Then, appealing to the results of Proposition 3, equations of motion 
of electromagnetic type with the desired property should be of the form: 

BV f)A a 

with 

A a = Tl{t)x b , and V = W(t,U T (t)x) + ±x T T 2 (t)x. (30) 

We illustrate this procedure by constructing a specific example for n = 3 by starting 
somehow in the middle of the converse procedure we just explained, i.e. in the represen- 
tation (|26|) . Consider the function 

W(t,y) = a{t)(yi -py 2 ) 3 , P constant. 
The general symmetric matrix S for which (|27p holds with this W, is given by 

(si k mp \ 

k 3i + (k/p)(l -p 2 ) m 
mp m S3 J 
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but to avoid expression swell, let us make the following choice for the five free constant 
parameters still at our disposal: si = S3 = k = p = 1, m = 0. Obviously the y-equations 
we are taking as starting point happen to exhibit a certain feature of decoupling already, 
they read 

y\ = -3a(t)(yi - y 2 f \ yi = 3a(t)(yi - y 2 ) 2 , m = 0. 
We choose U(t) to represent a simple rotation around de y\ axis: 

/I \ 

U(t) = cos 9(t) -sin0(t) 

\ sm9{t) cos0(i) / 

then r := —UU T becomes 

/ \ 
r= a 
\ — o / 

where we have put = All we need now to write down admissible x-equations 
of the form (|29p is the scalar potential V(t, x): it is given by 

V(t,x) = a{t){x\ — (cos^)x2 — (sin^)^) 3 — \a 2 [x^ + 

It is clear even without writing them down that the x-equations will be fully coupled. 
However, they have the alternative multiplier 

(1 cos 9 sin 9 \ 
cos 9 1 . 
sine 1 / 

The matrix g (just like S of course) has three distinct eigenvalues, namely 2, 0, 1. Hence, 
Proposition 2 proclaims that the x-equations should decouple completely in appropri- 
ate coordinates (and of course the partially coupled y-equations we started from just as 
well). It is a fairly straightforward exercise (at least with some computer algebra assis- 
tance) to compute a matrix P(t) which diagonalizes g{t) and to verify that the resulting 
transformation x = P(t)z does indeed have the predicted decoupling effect. But note 
with regard to the y-equations that the transition to the decoupling z-coordinates may 
have re-introduced velocity terms, resulting from the fact that the matrix P{t) which 
diagonalizes g need not be taken strictly orthogonal, and that P~ 1 U accordingly need 
not be constant. Of course, if we diagonalize the constant multiplier S of the y-equations 
via a constant matrix, we obtain another transition to a decoupled z-system without 
velocity terms. 
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